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Abstract 

We consider the effect of a partial coverage of quasar 
broad-line regions (QSO BLRs) by intervening H 2 - 
bearing clouds when a part of quasar (QSO) radiation 
passes by a cloud not taking part in formation of an 
absorption-line system in the QSO spectrum. That 
leads to modification of observable absorption line pro¬ 
files and consequently to a bias in physical parameters 
derived from standard absorption line analysis. In ap¬ 
plication to the H 2 absorption systems the effect has 
been revealed in the analysis of H 2 absorption system in 
the spectrum of Q 1 232+082 (see Ivanchik et al. 2010l 


Balashev et, al.ll201lli . We estimate a probability of the 


effect to be detected in QSO spectra. To do this we 
derive distribution of BLR sizes of high-z QSOs from 
Sloan Digital Sky Survey (SDSS) Data Release 9 (DR9) 
catalogue and assume different distributions of cloud 
sizes. We conclude that the low limit of the probability 
is about 11%. The latest researches shows that about a 
fifth of observed H 2 absorption systems can be partially 
covered. Accounting of the effect may allow to revise 
significantly physical parameters of interstellar clouds 
obtained by the spectral analysis. 
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1 Introduction 

An effect of the partial coverage occurs when an absorb¬ 
ing cloud does not fully cover an emission source. It is 
mostly probable at a situation in which an angular size 
of the source is comparable with a size of the absorbing 
object. The partial coverage is greatly important for 
spectroscopic studies including a choice of an appropri¬ 
ate model of the absorption-line systems (ALSs), which 
is necessary to derive physical parameters of an absorb¬ 
ing cloud from spectral analysis. In the simplest case 
of partial coverage the part of radiation from the emis¬ 
sion source passes by the absorbing cloud. It results 
in an additional residual flux at the bottom of absorp¬ 
tion lines registered in the spectrum with respect to the 
zero-flux level. This effect is more pronounced in the 
case of highly saturated lines. 

It should be emphasized that the most typical situ¬ 
ation for investigations of interstellar medium is when 
angular size of an absorbing cloud is much larger than 
an angular size of a background source which can be 
considered as a point-like object. In such cases the 
effects of the partial coverage are very unlikely and 
the observer uses simple absorption-line profiles defined 
by the optical depth r(A) according to the Bouguer- 
Lambert-Beer law. However, this approach ceases to 
be universal for studies of the ALSs in QSO spectra, 
especially for relatively compact H 2 clouds whose sizes 
can be comparable with sizes of QSO emission regions. 

QSOs are the powerful active galactic nuclei located 
at cosmological distances (up to z^7). Their high lu¬ 
minosity allows to detect their spectra and to study 
medium between QSOs and the observer via analysis 
of the ALSs imprinted in their spectra. Remind that 
the ALSs detected in QSO spectra can be classified in 
two general categories: intrinsic and intervening. In¬ 
trinsic absorption systems are related to a region in the 
close vicinity of the central QSO machine or the host 
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galaxy. Intervening absorption systems appear from in- 
tergalactic medium and interstellar medium of galaxies 
remote from the QSO. It was found that the effects of 
partial coverage, breaking the universal conception of 
a point-like source against an extended absorber, is a 
rather freque nt phenomenon for th e intrinsic ALSs in 
QSO spectra ( Petitiean et al. 1994 1. 

However, up to recently in the case of intervening 
ALSs the effects of partial coverage were believed to be 
unlikely. Actually, the projected sizes of typical inter¬ 
vening absorbers are much larger than quasar emission 
regions. The emission regions of QSOs in UV and op¬ 
tical range have been constrained by a size of < 1 pc. 
For example, reverberation mapping establishes the re¬ 
lationship between the BLR size and the luminosi ty and 


yields a BLR size of Rrt.r ~ 0.2 pc (|Kaspi et al 


2007 


C helouche a nd Daniel 20 12) for high redshift. luminous 
quasars. Differential microlensing all ows one to est i- 
mate the size of the BLR ~ 0.1 pc (jSluse e t al.1 12011 1. 
The observations of gamma-ray emission c onstrai n the 


size o f the jet extension to a few parsecs (jAbdo et al 
2010|) . On the other hand, majority of ALSs consti¬ 


tute the Ly-a forest, which are related to i ntergalacti c 
clouds with sizes well exceeding 10 kpc (Rauch! 1199811 . 
Sizes of diffuse atomic clouds inside intervening galax¬ 
ies are considered to imprint the damped Ly-a sys¬ 
tems (DLAs), or sub-DLAs, can be roughly estimated 
as 10 pc - 1 kpc. At last, the Lyman-limit systems 
(LLS), as well as Mg II and C IV absorbers, seems to 
correspond also to kpc-scale clouds in the halo of the 
intervening galaxies. 

Nevertheless, the partial coverage of Mg II system in 
the spectru m of the QSO, APM 08279+5255, was firstly 
reported by Ledoux et al. ( 1998ll . However, it had been 
found some later that this was very peculiar case of the 
partial coverage due t o a gravitational lensing of the 
backg round QSO (see Ledoux et al. 19981 : Ibata et al. 
19991 ). More definitely the effect of partial coverage 


has been revealed in the case of H I 21cm abso rptio n 
li nes measurements in the DLA system (Kanekar et al. 
2009). Actually, the emission regions of radio-loud 


QSOs at a wavelength 21 cm can well exceed a size 
of 10 kpc and typically consist of compact few pc cores 
and extended dilute envelopes. The effect of partial 
coverage essentially alters deri ved sp in tem peratu re T s 
of the neutral gas absorber ( Kanek ar et al. 20091) . 

Molecular bearing ALSs are another type of cur¬ 
rently observed absorption systems in QSO spectra 
which may have projected sizes comparable with a 
size of QSO emission region. These systems include 
molecular absorption lines (predominantly H 2 , in some 
cases - HD and CO) and refer to the diffuse and 
translucent molecular phases of the interstellar medium 


( Snow and McCall 200 ( 1 . Estimated number density 
of this system is > 10 cm~ 3 at the column density 
< 10 2n cm -2 , consequently these clouds constrained 
by pc or even sub-pc sizes (see Sect l3.ll) . 

The major part of QSO radiation is generated by an 
accretion disk (AD) around the black hole and forms a 
continuous spectrum. Characteristic size o f the AD _is 
considered to be < 3 x 10 -3 pc (e .g. iBlackburne et al 


201 It Jimenez-Vicente et al.J [2012). Thus this source 


may be treated as a point-like object fully covered by 
molecular clouds. On the other hand, there is still 
an extended region around the AD of relatively dense 
and warm partly ionized gas, called a broad-line region 
(BLR), where broad emission lines are formed, that 
leads to some additional enhancement of the contin¬ 
uum in QSO spectra. The angular size of this region at 
least an order of magnitude exceeds the typical angular 
sizes of ADs and may be comparable with sizes of H 2 
molecular clouds. Therefore an incomplete coverage of 
BLRs by molecular clouds situated on the line-of-sight 
between the QSO and observer becomes quite probable. 

The partial coverage of a BLR by t he intervening 
H 2 absorption cloud was reported by Iwmchil^etjil 


( 201C ) and investigated in details by Balashev et al 


(2011). It was demonstrated that the H 2 -bearing cloud 
( z a bs = 2.34) covers completely the QSO 1232+082 
(z em = 2.57) intrinsic continuum source, but only a 
part of the BLR. 

There are also other possibilities for the partial 
coverage effect: (i) jet emission can be scattered/re¬ 
emitted by medium in the quasar host galaxy, which 
is remote from the AD and BLR. In this case the 
observer detects the partial coverage of a continu¬ 
ous spectrum. This is the most probable explana¬ 
tion of the recently revealed partial coverage of the 
Q 0528-250 conti nuous spectrum by H 2 absorption sys¬ 
tem at z=2.811 ( Klimenko et al. 2015 ). (ii) The cal¬ 
culation of galaxies evolution shows that a galaxy 
at high redshifts can have several supermassive black 
holes (SMBHs; binary, triple or quadru ple systems 


noies (bivl-tsils; binary, triple 01 +quadruple systems 
of SMBHs, e.g. lKulkarni and Loebl 2012 ). in partic¬ 


ular, one of SMBHs can remain unscreened by an ab¬ 
sorbing molecular cloud, (iii) The observer can reg¬ 
ister UV emission of a QSO host galaxy much more 
extended than a BLR. (iv) The partial coverage can 
be caused by gravitational lensing (e.g. the case of 
APM 08279+5255) of a galaxy or a galaxy cluster lo¬ 
cated between a quasar and intervening clouds. 

The effect of partial coverage leads to serious changes 
of parameters derived from an analysis of Q SO spec- 
tra (e.g. column densities of absorbers, llvanchik et al 


201(1 Balashe v et al. feoill ; iKlimenko et alj |2015|j . So it 

is worthwhile to estimate a probability to find the par- 
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tial coverage of a QSO BLR by an intervening H 2 ab¬ 
sorption system. Other appearances of the partial cov¬ 
erage (listed above), apparently, lead only to increase 
the probability. 

Following lBalasliev et al. ( 2011 ). one can introduce a 
flux coverage factor f as the ratio of a light flux compris¬ 
ing the absorption systems imprinted in its spectrum to 
the total flux from a quasar. In the present paper we 
calculate a probability to reveal the partial coverage of 
a QSO BLR located at some redshift z q by H 2 ALS at 
redshift z c (cloud) with coverage factor / < f 0 (/o is a 
fixed value, see below). We consider an arbitrary angu¬ 
lar distance between line-of-sights to the centres of the 
cloud and BLR, as well as an arbitrary ratio n of the 
cloud to BLR transverse sizes. 

The paper significa ntly develops our previous work 
(Ofengeim et al. 120131 ) and is organized as following. 
At first, we proceed from an oversimplified assumption 
of uniformity of the radiation flux from the whole BLR 
region ignoring the luminous point-like AD in the cen¬ 
tre of the QSO. In Section[2]we introduce a geometrical 
coverage factor g = / and consider the main parame¬ 
ters controlling g. In Section [3711 we calculate a prob¬ 
ability to reveal factor g < go-, where go is also a fixed 
value (an analogue of /o), for certain distributions of 
the main parameters. In Section 13.21 we discuss simi¬ 
lar probability but for more realistic coverage factor / 
determined as the ratio of radiation fluxes, taking into 
account nonuniformity of the radiation flux and includ¬ 
ing the intense AD radiation. Sections EHfl andHloresent 
a special consideration of the dependence of the factor 
/ on the redshifts z q and z c . In Section [5] we discuss 
obtained results and their relation to observations. 


2 Geometry of partial coverage 

In this work we use a toy model to describe the emission 
regions of QSOs and absorbing clouds. The geometry 
of this model is presented in Fig. 11121 The transverse 
sizes of the BLR and the cloud are designated as l q 
and l c , the angular size of the BLR is designated by 
0 q (hereinafter angles are in radians), a correspond¬ 
ing angle for the cloud can be denoted as 9 C (is not 
shown in Fig. [I]). The light from the BLR propagates 
through the Universe and it is registered by the ob¬ 
server. The observer detects an emission within a light 
cone restricted by an angle 9 q or a solid angle 9 g com¬ 
prising the whole flux from the BLR. Some part of light 
from the BLR is partly shielded by the cloud in such 
a way that molecules inside the cloud imprint a set of 
absorption lines (or ALS) in the initial spectrum. This 
part is bounded by a solid angle (of coverage) 9 co „. 


The rest light can pass by the cloud and come to the 
observer without formation of ALSs in its spectrum. 
This flux of radiation may be approximately estimated 
as a flux comprised by a solid angle Q unC ov = 9 g — 9 co „. 
In result the observer detects a complex radiation from 
the QSO with spectrum integrated over angles. 

Let us determine a geometrical coverage factor g as 
a simplest way to estimate factor /. If we assume that 
fluxes from covered and uncovered parts of the BLR 
are uniformly distributed within their solid angles, then 
factor g may be estimated as a ratio of the solid angle 
fl C ov to the whole solid angle Q q : 


fL 


9 = 


9„ 


( 1 ) 


where at 9 q -C 1 we can write f l q = n9^. We imply 
that g = f in the case of pure geometrical estimations. 
Actually, the factor g may essentially differ from / be¬ 
cause it does not take into account nonuniformity of 
fluxes within the solid angles (see Section [372]) . 

Thereafter it is convenient to express all angular val¬ 
ues in units of 9 q and all solid angles - in units of f l q . 
Thus a relative angle between directions to the centres 
of QSO and cloud (angular deviation) is determined as 



( 2 ) 


and the relative angle size of a cloud becomes 


P = 


9c 

V 


( 3 ) 


An angle size 9 of an object with a transverse 
size l at the cosmological redshift z may be d eter¬ 
mined by a st andard way (e.g. IZeldovich and Novikov : 


Kavser et alj 19971) as 


9 = 


l 


2 D a {z)' 


( 4 ) 


where D A {z) is an angular diameter distance. In the 
standard ACDM cosmological model the value D A (z) 
to a cosmolog ically distant object at the redshift z is 
(e.g. IZeldovich and Novikovt iKa vser et al, 1997h 


D a (z) = 


Hr 


— t 

1 + z Jo 


dz' 


0 \/9 m ( 1 + z')3 4 - 9a 


( 5 ) 


where c is the speed of light, H 0 = 100 h km s -1 is the 
present Hubble constant, c/H 0 ~ 2.998 h~ x Gpc; 

9 m is the dimensionless matter density parameter 
and 9 a = Ac 2 /( 87 tG) is determined by the cosmo¬ 
logical constant A, G is the gravitational constant; 
hereafter we assume that 9 m ~ 0.3 and 9 a = 0.7 
1 Planck Collaboration et ah 2014) 
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Fig. 1 Schematic illustration of an absorption cloud at a transverse size l c (dark grey circle) situated between the observer 
and QSO BLR with a transverse size l q (white ellipse in the centre of QSO); the wider ellipse around BLR (light grey) 
symbolizes a QSO host galaxy; 6 q is an angular size of BLR, AO is an angle between two line-of-sights from the observer to 
the centres of QSO and cloud; z q and z c are redshifts of the QSO and cloud; 9l q is a solid angle (light cone) of the whole 
BLR radiation flux, Q cov is a solid angle of a part of BLR flux passed through the cloud. Radiation from the QSO without 
traces of hQ-absorption systems in its spectrum is shown as a white area of the light cone and radiation containing the 
absorption systems - as a light grey area. In general case the light cone is curved due to the expansion of the Universe. 


Using Eq. (j4]) one can rewrite Eq. m as 


p{K, z q ,z c ) 


Da(z q ) 
Da(z c ) ’ 


( 6 ) 


where the relative angular size of the cloud is expressed 
through the ratio of transverse sizes of the cloud and 
QSO 



( 7 ) 


The coverage factor g is a function of p — p(n; z q , z c ) 
and 5. Fig. [2] demonstrates four basic types of relative 
positions of two circles conventionally referred to the 
BLR and cloud projected on the sky. 

These types corresponding to special dependence 
g = can be formalized as follows: 

(i) Full coverage: 1 + 5 < p. All flux from the QSO BLR 
goes through the cloud, g = 1. Up to recently, only this 
case has been considered in literature in application to 
intervening LU clouds. 

(ii) Partial ring-shaped coverage'. p+S < 1. The angular 
size of the BLR exceeds the angular size of the cloud 
plus the angular deviation. In this case g = p 2 . 


(iii) Crescent-like coverage : \p — 1| < <5 < p + 1. After 
some calculations one can obtain 


, n ^ P 2 P 2 +S 2 - 1 

g(p,d) = s{p , d) = — arccos-—- 

7T ZpO 


1 

H— arccos 

7T 



1 + S 2 — p 2 
2 6 

-(6-ims+iy 


( 8 ) 


(iv) Full noncoverage: 1 + p < S. Relative angular 
distance is large and the cloud remains unobservable in 
the QSO spectrum. In this case 5 = 0. 

One can summarize all cases (i)-(iv): 


g{n,S; z q ,z c ) 


i; « 

P 2 {k\ z q ,z c )\ (ii) 

s[p(n; z q ,z c ),8\; (iii) 
0 ; (iv) 


(9) 


3 Distribution of coverage factor 

3.1 Geometrical coverage factor 

Let us calculate the cumulative distribution function 
(CDF) P{g < g 0 ), which describes the probability to 
detect an absorption system in QSO spectrum with ge¬ 
ometrical coverage factor g not exceeding a value go. 
Since one can reveal the cloud only detecting absorp¬ 
tion lines in QSO spectra, the case of full noncoverage 
g = 0 is equivalent to the absence of an observed cloud. 
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5+l<p p+8<l |p-l|<5<p+l p+l<5 


Fig. 2 Four types of coverage (grey regions) of QSO BLRs (white circles) by clouds (black circles) in the observer reference 
frame: (i) full coverage, <7 = 1 (widely accepted case, white circle coincides with grey one); (ii) partial ring-shaped coverage, 
0 < g < 1 (black circle coincides with grey one); (iii) crescent-like coverage, 0 < g < 1; (iv) full noncoverage , g = 0. All 
angular sizes are represented in units of 9 q \ the angular size (radius) of the BLR is 1, the radius of the cloud is p, and the 
angular deviation is 5 (see text). 


Therefore we can put (rather conventionally) that the 
probability of full noncoverage is P{g = 0) = 0. 

More specifically, to detect an absorption system in 
a spectrum it is necessary that a cloud totally cover 
the QSO AL0. In contrast, if the cloud partially cover 
only BLR and does not cover AD, tQ molecular ab¬ 
sorption lines associated with the cloud in the spec¬ 
trum will be detectable only in surroundings of the top 
of quasar emission lines blueshifted of Lya (e.g. C ill 
or Ly/3+Ovi, see Fig. 0. However, in such a case the 
effect of the partial covering would be difficult to iden¬ 
tify, since only several weak absorption lines of the H 2 
ALS could be expected. These lines would be hardly 
detectable against the numerous Ly-alpha forest lines 
in the spectrum. It gives an additional restriction on 
the angular deviation: 5 < p, or P(6 > p) = 0. 

Let us add, that the effect of the partial coverage 
can be observed also for Civ emission line (redshifted 
of Lya) using C 1 (neutral carbon) absorption lines as¬ 
socia ted with H 2 ( Ivanchik et al. 20ld Balashev et al.l 
201ll) . but this case is not considered in this paper. 


The CDF P(g < g 0 ) can be determined through the 
distributions of two random values n' and 5'. Let us 
determine the probability density functions (PDFs) for 
these two distributions: p K (re) andp 5 (d| re), respectively, 
the latter value being the conditional PDF of the occur¬ 
rence of S' = <5 at a given value re/ = re. Following to the 
standard definitions of the Probability theory one can 
introduce the relations between these PDFs and differ¬ 
ential probabilities: P(re < re' < re+dre) = p K (re) dre and 
P(S < S' < 5 + d<5) = ps(6\ re) d<5, and as a consequence 


1 Remind that the AD generates the most part of continuum ra- 
diation 


to write 


P{g < go) = 

/•OO pCX 

Jo Jo 


pa(6| k)p k (k) Q[go - g(re,<5)] d£dre. (10) 


where 0(a;) is the Heaviside function. One can optimize 
calculations using the reciprocal relation re(< 7 , 6) given 
in Eqs. (l26l) and (l27l) of Appendix. Then integration 
over re in Eq. dll should goes from 0 to n(go, S). 

For simplicity, we assume the uniform mutual distri¬ 
bution of QSO and cloud centres on the sky and get 
ps(5) oc 6. To obtain a normalization function for ps(S) 
we use the condition S < p of the reliable detection of 
H 2 ALS (see above). In result we obtain 


Ps(<5| re) = 


26 


[p(re)] ; 


0[p(re) - <5], 


( 11 ) 


It was chosen that z q = 2.57, z c = 2.33 in all cal¬ 
culations employing Eq. d 6 j) for p(re, z q , z c ), he. they 
refer to the case of th e partial coverage of Q 1232+082 
(iBalashev e t al.l 12011 ). The dependence of results on 
the redshifts will be investigated in Section 13.31 

The PDF p K with respect to re may be represented 
as 

pOO 

Pk.{k) = 2 / Pdoud(2.Kr) Pqso(r ) rdr, (12) 

Jo 

where p c ioud and p qso are PDFs for linear sizes of clouds 
l c = 2rer and radii r = l q /2 of BLRs, respectively; in 
derivation of Eq. m it is used that d l c = 2 rdre; factors 
dre are reduced in both sides of Eq. dm 

To obtain p n a n(r) w e exploit an empirical equation 
by Kaspi et ah ([2005 ) between the radius r and QSO 
luminosity 


10 light-days 


= AL ± 


(13) 
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Fig. 3 Left panel: Histogram (hatched bars) of QSO-luminosities sampled at a wavelength of 1450 A from the data of 
SDSS DR9, QSOs are chosen under condition of the presence of DLAs in their spectra; right panel: probability density 
function (PDF) of BLR radii, hatched histogram bars form a distribution obtained with using Eq. (1131) and the data shown 
in the left panel; solid line is the resulting lognormal approximation. 



Fig. 4 Left panel: Three types of cloud sizes, l c , distribution. Short-dashed curve is the uniform distribution within an 
interval of 0.2 < l c < 15 pc, solid curve is the lognormal distribution on a condition that the values l c with 90% probability 
belong to the same interval 0.2 < l c < 15 pc, long-dashed curve is the lognormal distribution on a condition that the most 
probable value is l c = 0.2 pc (see text). Right panel: the cumulative distribution function P(g < go) — the probability 
to detect absorption systems in a QSO spectrum (represented at z q = 2.57 and z c = 2.33) with the geometrical coverage 
factor g < go- The upper argument, fo (flux coverage factor), determined by a scaling fo = 0.3</o + 0.7 in such a way that 
the cumulative distribution functions P(f < fo) remains unalterable. Types of lines are the same as in the left panel. 
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where L = ALa/(10 44 [erg/s]) is the normalized QSO 
luminosity at a fixed QSO-restframe wavelength A, A 
and B - numerical coefficients. To estimate QSO lu¬ 
minosity di stribution we used QS O spectra from SDSS 
DR 9, (see Ahn and et ajJ 120121) . It is commonly ac¬ 
cepted that due to high column density H 2 -cloud asso¬ 
ciated with large amount of neutral hydrogen, i.e. asso¬ 
ciated with D LAs. We used the samp le of DLAs bear¬ 
ing QSOs ( Noterdaeme and et al.ll2012 ) from the whole 
SDSS DR9 QSOs sample. Choosing in Eq. (fl3l) lumi¬ 
nosities at the QSO-restframe wavelength A = 1450 A 
and quantities A = 2.12, B = 0.496 (errors omitted) 
we can calculate a required distribution p qso (r). 

Actually, Fig. [3] shows the distributions of QSO 
(DLA) luminosities at A = 1450 A (left panel) and 
Pqso(j* ) calculated for BLR radii (right panel). One can 
see that the lognormal approximation 




Pqso{r) = 


\f2irv 


(14) 


where w = 0.426 ±0.003 and ro = (0.1442 ± 0.0004) pc 
obtained by the least square method, proves to be good. 

The situation with a distribution Pdoudih ) for H 2 
molecular clouds is still more uncertain. Absorption 
systems of H 2 detected at high redshifts are referred to 
the cold neutral medium, which has a lower limit ofthe 
number density about 1-710 cm -3 (e.g. see Lisztll2002 ). 
It gives l c < 3 -7 30 pc for typical measurable column 


analysis of I D absorption system in the spectrum of 
Q1232±082 ( Balashev et al. 2011 1 yields l c ~ 0.2 pc. 
Therefore we employ for our estimation two modeled 
distributions: (1) the uniform distribution within an 
interval [0.2, 15.0] pc 


^L(A) = r i -B(A-a)B(^A), 


b — a 


(15) 


where a = 0.2 pc and b = 15.0 pc. One can expect that 
this model gives overestimated sizes of the clouds and 
consequently underestimated values of P(g > go)- (2) 
The lognormal distribution 

(In l c 

(in 7p ) 

P^loud^c) = -> ( 16 ) 

where parameters u and lo are choosing also for two 
cases. The first one is u = 1.2 and Iq = 0.8 pc, it 
corresponds to the fixed maximum of the distribution 
Eq. (fTTH) or to the most probable size at l c = 0.2. The 
second one satisfies to the condition that 90% clouds 
fall into the interval [0.2,15.0] pc. In that case u = 1.3 
and Iq = 2.0 pc. 


The left panel in Fig. [4] demonstrates three types of 
the PDF for the cloud sizes l c defined by Eqs. (fl5l) and 
(11611 with two sets of parameters u and Iq. The right 
panel shows the integral probability P(g < g 0 ) as func¬ 
tions of go calculated with the use of Eqs. m,m, and 
m at all three types of distribution p c ioud(lc = 2 nr) 
discussed above. Let us emphasize the strong depen¬ 
dence of the probability distribution on the distribu¬ 
tion of cloud sizes. Nevertheless, even the most conser¬ 
vative estimation of the probability to reveal g < 0.98 
proves to be « 15%, where go = 0.98 corresponds to the 
minimal measurable n oncoverage factor: 1 — go A 0.02 
( Klimenko et al. 2015 1. 


3.2 Flux coverage factor 


In contrast with the oversimplified geometrical coverage 
factor g one can introduce the factor / as a ratio of 
radiation fluxes, which is more relevant to a spectral 
analyses: 


n _ *COV 

J — J ! 

J-qso 


(17) 


where I cov is the flux of the quasar radiation passed 
through the cloud, I qso is the total flux. This value 
may be called as a flux coverage factor. In the case of 
uniform distribution of the flux over a solid angle Q q 
we have I cov oc D co „, I qso oc and the flux factor / 
would be equal to g. 

Let us regard, some arbitrary, that the flux from the 
AD characterized by a continuum spectrum constitutes 
about 70% of the whole intensity. The rest 30% of the 
flux originates from the BLR (Fig. [5j and represents a 
set of broad emission lines which also raises a level of 
continuum in some parts of the spectrund 

If the AD is covered totally but the BLR is covered 
partially (left panel in Fig. [5]), the partial coverage effect 
(/ < 1) will be detected only in those absorption lines 
which are located around the top of QSO emission lines. 
Outside the wide emission lines the partial coverage is 
not detectable, that is equivalent to / = 1. Within 
the regions of emission lines the full covering of the AD 
flux yields only a part of the covering factor / = 0.7 
and the rest part of / is provided by the geometrical 
coverage factor g related solely to the BLR. As a result 
we assume that the flux coverage factor / is determined 
by a simple equation: 


/ = 0.7+ 0.3$. 


(18) 


2 The fact that a continuum source can be fully associated with 
the AD is proved b y observations of rapid continuum variability 
(see Petersonlll993l and references therein) 























Fig. 5 Model representation of radiation-flux nonuniformity. Left panel: the accretion disk of QSO is covered by a cloud 
totally, Ad < 9 cloud. (5 < p); right panel: the accretion disk of QSO is uncovered by a cloud, A 9 > 9 c i ou d (5 > p). The 
flux passing by the cloud contains the whole flux of the accretion disk and more than a half of flux emitted by the BLR. 
In both case the BLR is covered only partially, but the case represented on the right panel is unlikely to be detected (see 
Section lTlI for details). 


Then, the calculation of the CDF P(f < fo) is re¬ 
duced to calculation of P(g < (fo — 0.7)/0.3), i.e. to 
pure geometrical approach. A scaling /o = 0.7 + 0.3go 
applied to the upper horizontal axis in the right panel 
of Fig. [Hallows one to transform the dependencies on g 
— P(g < go) into the dependencies on / — P(f < f Q ). 

The most conservative estimation of the probability 
to find a partial coverage at the minimal level of the 
noncoverage factor 1 — / 0 > 0.02 or f 0 < 0.98 (e.g. 
Klimenko et al. 120151) yields 11% (lower curve in the 
right panel of Fig. a. Note, that to detect the residual 
flux at the bottom of lines (which typically include sev¬ 
eral pixels) at the level ~ 2 % it is necessary to obtain 
high resolution spectra with high signal to noise ratio, 
S/N ~ 30 — 50. Spectra registered with less S/N allow 
to detect only higher minimal level, e.g. 1 — /o > 0.1 
or fo ^ 0-9- It gives the lowest level of the probabil¬ 
ity: ~ 5%. However, we suppose that the middle solid 
curve in Fig. [4] is more realistic. It gives the probabil¬ 
ity estimation of the order of 15% to reveal the partial 
coverage with / < 0.9. Thus, the effect of partial cov¬ 
erage of QSOs by intervening H 2 ALSs characterized 
by noticeable probability to be revealed at the spectra 
analyses. 


3.3 Dependencies of /-distribution on redshifts 


0.6 
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Fig. 6 The cumulative distribution P(fo', z q ) = P(f < 
.fo',z q ) (see text) as a function of QSO redshift, z q , when 
a cloud redshift z c = z c (z q ) is maximally remote from z q 
(see Eq. ©). Types of lines are the same as in left panel 
of Fig. [T] Each pair of the curves corresponds to the same 
distribution of l c m , the higher curve of a pair is plotted at 
fo = 0.98, while the lower one - at fo = 0.90 (numbers near 
the curves). Crossed points correspond to the special case 
of P(fo\ z q ) calculated for Q1232+082 at z q = 2.57 and the 
H 2 -cloud at z c = 2.33 registered in its spectrum. 



Let us consider dependencies of the probability distri¬ 
bution on the redshifts z q and z c . Taking into account 
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dependence of the distribution p^(<5|ac) on z q and z c , de¬ 
fined by Eqs. (HD and ©, one can modify Eg. [TUland 
obtain CDF as 

P(f < fo;z q ,z c ) = 

roo rK,[(fo—0.7)/0.3,S,z q ,z c ] 

/ / Ps(S\k, Z q , Z^p^K.) d8 dK, 

Jo Jo 

(19) 

where according to Eq. m fo varies between 0.7 and 
1.0, and the function k is given in Eqs. m and USD. 

Firstly, the dependence of the CDF on redshifts are 
expected to be small, since the angular diameter dis¬ 
tance (see Eq. ([5])) weakly depends on the redshifts 
within an interval z c , z q ~ 2 — 4. Let us consider 
an H 2 -cloud, which is remote from the QSO by the 
maximal distance consistent with the very possibility 
to detect a set of H 2 absorption lines in surrounding of 
a broad emission line and consequently to detect the 
partial coverage. 

One can roughly formulate this condition of max¬ 
imal removing as (1 + 2 C )1094A > (1 + z q ) 977 A, 
where 1094 A is wavelength of Ll-0 P(l) transition and 
977.02 A corresponds to the C ill prominent emission 
line in the restframe (see Fig. 0, since most of satu¬ 
rated H 2 absorption lines are blueward in relation to 
1094 A and all prominent emission lines of QSO spec¬ 
tra are redward in relation to 977.02 A. Note that we 
conservatively take Ll-0 P(l) line as a boundary, since 
the most long-wavelength H 2 transitions correspond to 
LO-O band at A ~ 1108— 1110 A (see Fig.0 frequently 
appears as a unsaturated absorption lines and can not 
be used for confident detection of the partial coverage. 
Therefore, to obtain partial coverage of the QSO BLR 
by a H 2 -cloud one needs that redshifts of H 2 ALSs obey 
the condition 

0.89 z q — 0.11 < z c < z q , (20) 

where we put 977/1094 ~ 0.89. The right part of con¬ 
dition just indicates that H 2 ALS is located between 
QSO and observer. 

Fig. [^demonstrates dependencies P(/o; z q ) = P(f < 
fo; z q , z c (zq )) on z q at two fixed fo = 0.98 and 0.90 and 
z c (zq ) determined by minimal values of Eq. (l20l) . One 
can see that the probability to reveal the flux coverage 
factor 0.7 < / < /o weakly depends on z q . More¬ 
over, the probability P(fo, z q , z c ) weakly depends also 
on z c . Actually, it is demonstrated by crossed points 
in Fig. [G] calculated for the QSO 1232+082 (z q = 2.57) 
whose spectrum contains the H 2 absorption system at 
z c = 2.33. The crossed points lie closely to the corre¬ 
sponding curves P(/o; z q ), while the lower boundary z c 
at z q = 2.57 according to Eq. (l20l) is z c > 2.17, i.e. the 


difference does not affect noticeably the value P(fo;z q ). 
Thus we can disregard approximately the dependence 
P(fo; Zq) on z c . 

4 Notes on observations of the partial coverage 

effect in QSO spectra 

As an illustration Fig. [7] displays typical spectrum of the 
QSO in the restframe. Specifically we used the quasar 
spectrum J 0003-2323 at z q = 2.28. It was observed 
under program 166.A-0106 (PI: Jacqueline Bergeron). 
The data are in free access from the ESO database!. 
One can see strong Ly a emission line, a number of ab¬ 
sorption lines of the Lyman-alpha forest imprinted in 
a continuum, damped Lyman-alpha absorption (DLA) 
system at Adla = 1215 A(1 + z c )/(l + z q ), broad emis¬ 
sion lines OlV, Ly^ and Cm, and schematically a set 
of H 2 -molecular absorption bands. We regard (see Sec¬ 
tion GED that H 2 -bearing clouds are supposed to be 
associated with DLA systems. 

Fig. [8] shows the location of known H 2 absorption 
systems at high redshifts on the z q — z c plane. All H 2 
systems capable to cause the partial coverage of QSO 
BLRs can be placed mainly between two linear bounds 
determined by Eq. m- Fig. [5] represents also prelimi¬ 
nary results of a systematic search of H 2 -systems with 
traces of the partial coverage effect (Klimenko et al, in 
prep.). All five candidates marked by crosses in Fig. [8] 
have to be verified in further analysis. 

5 Conclusions 

In this paper we estimate the probability of an incom¬ 
plete coverage of the QSO BLR by the absorbing H 2 
molecular cloud. The key value of our analysis is a cov¬ 
erage factor / - a ratio of the flux of radiation which 
goes through a cloud to the total flux could have reg¬ 
istered from the QSO without absorption by the cloud. 
We calculate the probability to obtain different cover¬ 
age factors at given cosmological redshifts z q and z c 
of the QSO and cloud, respectively, employing the dis¬ 
tribution of BLR sizes based on the data of the SDSS 
DR9, and a few model distributions of cloud sizes. 

According to Fig.[4]the most conservative assessment 
of the probability to reveal the factor / < 0.98 (obser¬ 
vational minimum of noncoverage factor 1 — / > 0.02) 
is 11%. Thus one can expect at least three cases with 
/ < 0.98 from 27 H 2 absorption systems detected in 
QSO spectra up to now. The preliminary results of 


3 htt p: / / archive .eso.org 
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CIII 977A LyP +Ovi 1025A Lya emission 1215. 67A 
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Restframe wavelength, A 

Fig. 7 Spectrum of quasar J 0003-2323 at z q = 2.28 shown in its own restframe. There are marked Ly a emission line, 
DLA system at Adla = 1215A(1 + z c )/( 1 + z q ), broad emission lines Oiv, Ly^ and Cm, Grey vertical bands indicate 
the wavelength regions in surroundings of the broad emission lines where the effect of partial coverage could be mostly 
pronounced. In the upper left part of the figure it is shown also the series of typical absorption bands of the H 2 molecules. 
The spectrum is not properly flux calibrated, since it was obtained using Ultra Violet Echelle Spectrograph at VLT. For 
data see the footnote in the Section [4] 

the systematic search of the partial coverage effect in 
quasar spectra with imprinted intervening H 2 absorp¬ 
tion systems give five cases, that is consistent with the 
prediction. 

Systematic search for the partial coverage effect in 
the H 2 absorption systems will allow to obtain restric¬ 
tions on characteristic scales of the H 2 clouds. For in¬ 
stance, relatively high frequency of revealing the effects 
of incomplete coverage with 0.7 < / < 1.0 would mean 
that cloud sizes are comparable with the BLR sizes, 
i.e. ~ 0.2 — 0.3 pc. Relatively low frequency of such ef¬ 
fects in QSO spectra would mean that the typical cloud 
size essentially exceeds the BLR size as it was widely 
accepted previously. 

In this paper we consider a simple model, where ab¬ 
sorption and emission regions are represented by spher¬ 
ical shape clouds. In case of asymmetrical and/or 
patchy structure of the absorption as well as emission 
regions the situation can be sufficiently complicated by 
increasing the number of parameters and uncertainties. 
However, it seems that estimate of the probability of re¬ 
vealing partial coverage - the main result of this paper, 
will not significantly changed. Our understanding of 
the structure, distribution of sizes and physical condi¬ 
tions of cold interstellar medium at high redshift would 
be constrained by the further observations and exten¬ 
sion of the statistics of the BLRs partial coverage by 
H 2 absorption clouds. 

Acknowledgements This work is supported by the 
Russian Science Foundation, grant 14-12-00955. 
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Fig. 8 Known high-redshifted H 2 absorption systems on 
the z q — z c plane. Empty circles indicate absorption sys¬ 
tems, where partial coverage has not been detected up to 
now (22 objects); crosses mark systems for which a prelim¬ 
inary analysis evidences in favour of an occurrence of par¬ 
tially covered BLRs (5 objects) in QSOs spectra. The upper 
solid line corresponds to z c = z q and determines the highest 
positions of H 2 systems on the z q — z c plane; the dashed line 
is z c = 0.89 z q — 0.11 (see Eq. (1201) 1 and gives the lowest z c at 
fixed z q when partial coverage by i nterveni ng cloud can be 
detected. The data was taken from llvanchik et alj (l2015h . 
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In Appendix we obtain an analytic approximations of 
g{n,5) determined in Eqs. © and © and represented 
as a set of sections continuously sewed together in sur¬ 
roundings of definite boundary lines k(S) (see below) 
at fixed z c and z q . Such an approach allows us to de¬ 
rive an appropriate set of simple reciprocal functions 
n(g, S) and to use them in the integrations of Eqs. (fTUl) 
and (fl9l) . 

Firstly, using Eq. © one can define two basic bound¬ 
ary functions on S at the plane S — n: 


K! 


D(z c ) 

D(z q ) 




— tJ-4(<5 + 1)- (21) 

D(z q ) 


These functions correspond to the boundaries p = |<5—1| 
and p = <5 + 1, respectively, at the plane S — p, which 
separate the crescent-like coverage region (iii) from the 
regions (ii), (iv) and (i) (see Fig. 0. 

Then, to obtain simple approximations within the 
region (iii), we introduce two additional interfacing de¬ 
pendences on 5 

2 1 12 

k 1/3 — g K l + g K 2 K 2/3 — -Ki + -K2, (22) 


and using Eqs. © and © calculate corresponding ac¬ 
curate functions 


Si/3 = s(« 1/3, S) 32/3 = 3(« 2/3, <5)- (23) 

Note, that all defined quantities depend on z q and z c , 
which are incorporated in the dependences k(S). 

To obtain the domain-like approximations for g{n, S) 
it is convenient to introduce two basic zones: (a) with 
S > 1 and (b) with <5 < 1. 

In the zone (a) at 6 > 1 we have: 


g{n,S) = 


0; 

0 < K < K\ 

/ X 3/2 

si/a ur, 3-0 ; 

K\ < K < K 1/3 

3l/3 + 

(32/3 3i/ 3 ) k 2/3 -k 1/3 , 

Kl/3 < K < ^2/3 

1 (1 32/ 3 ) (« 3 -« 2/3 ) > 

«2/3 < K < H 2 

i; 

K 2 < K, 


This manuscript was prepared with the A AS IAT^K macros v5.2. 


where five sections consequently correspond (from top 
to bottom): 1-st section to the region (iv) at p < S — 1; 
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2-nd, 3-d and 4-th sections to (iii) at <5 — 1 < p < <5 + 1; 
5-th section to (i) at <5 + 1 < p. 

In the zone (b) at 5 < 1 we have: 


g(n,6) = 
p 2 (k); 


0 < K < K\ 


9i/a~ A _ 3/2 «?/3 A - k i 2 9i/3 

3/2 ^3/2 K ' 3/2 3/2 ’ 


K l/3 K 1 


3l/3 + (S2/3 “ 3l/3) 


K l/3 K 1 

Kl < K < K1/3 

k—k, 1 / 3 , 

«2/3-«i/3 ; Kl /3 < K < «2/3 

3/2 


i; 


/ V J/Z 

1 - (1 - 32/3 ) (k^ K2 K /3 J 5 «2/3 < K < «2 

K 2 < K, 


(25) 


where A = (1 ^ <5) 2 , five sections consequently corre¬ 
spond: 1-st upper section to the region (ii) at p < 1 — S', 
2-nd, 3-d and 4-th sections to (iii) at 1 — <5 < p< l-l-<5; 
5-th section to (i) at 1 + S < p. Note that a differ¬ 
ence between the exact coverage factor g(n,5) and the 
approximations given by Eqs. (l24l) and (E3 does not 
exceed 0.02. 

Finally, we can obtain an approximation for the re¬ 
ciprocal function K,(g,S), treating the same two zones: 
(a) with <5 > 1 and (b) with S < 1. 

Thus in the zone (a) at S > 1 we have: 


In the zone (b) at 5 < 1 employing Eq. d6]) we have: 

K(g,S) = 

' 0 <g< a 

k i/3 (3 - A) + k 3/2 (si/3 -g) 

3i/3 — A 

‘ A < g < g 1/3 

*4/3 + ( K 2/3 - k 1/3) g 9 2 / 3 -gl /3 ; 3l/3 — 3 < 32/3 

( 1- \ 2/3 

K 2 - (k 2 - K2/3) [j^j;) ; 32/3 < 3 < 1, 

(27) 

where as in Eq. (l25l) A = (1 — <5) 2 , and the first upper 
section corresponds to the region (ii), all three others 
to (iii). 


n{g,6) = 


+ («i/3 - ki) ( 5*7) 


2/3 


Kl/3 + (K2/3 — K1/3) 


g-g 1/3 . 
92/3—31/3 ’ 

2/3 


0 < 3 < 3i/ 3 
3l/3 < 3 < 32/3 


„ K2 - (k 2 - K2/3) (137^7 J ; 32/3 < 3 < 1, 


(26) 


where all three sections correspond to the region (iii) 
and its boundaries with the regions (iv) at g = 0 and 
(i) at g = 1. 


















